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MOMENCLATURE 


Specific heat of steel 

Diameter of ingot 

Heat of fusion of steel 

Surface coefficient of heat transfer, ingot to mold 
Surface coefficient of heat transfer, water side 
Mold wall thickness 

Mold length 

Heat transferred, Btu/hr.,ft@ 

Heat transferred, Btu/hr. 

radius of ingot 

radius of mold, water side 

average surface temperature of ingot within mold 
average water temperature in mold 

mean temperature of ingot cross-section at mold exit 
pouring temperature 

Composite heat transfer coefficient 

velocity of ingot 

thermal diffusivity 

density of steel 

Nusselt number, water 

Reynold's number, water 


Prandtl number, water 
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Analysis of contimous casting involves many phases, two of which 
are described in this thesis. The first investigation consists of 
determining the temperature distribution in a moving ingot of metal 
which is solidified and cooled sufficiently in a water cooled mold 
to permit the contimous processing of the resulting billet. in 
the second phase there has been established a relationship for maximum 
water flow rates for heat transfer and their association with cor- 
responding coefficients of heat transfer. 

The analysis has shown that the rate of heat transfer is primarily 
a function of the surface coefficient of heat transfer between the 
ingot and the mold. To increase this rate some means mst be provided 
to increase this coefficient. Addition of oil along the boundary has 
been suggested which would not only increase this surface coefficient 
but will also decrease the frictional drag, 

The results, although only approximate, are provided as a guide 


for future work. 
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INTRODUCTION 


Generally, the analytic mode of approach to the solution of a 
problem also requires experimentation. However, lack of facilities and 
time limitations precluded experimentation for verification of any 
results in this thesis. 

The primary endeavor in this investigation is to assist the 
experimenter in his work by forming limits of practicability and 
providing mathmatical solutions to important phases of the problen. 
Results are only approximations, consistent with assumptions necessary 
to perform the operations, but they can be used as a guide for future 
work. 

Treatment of the entire problem of contimous casting is very 
long and complex, and to the best knowledge of the writer no work has 
been done analytically along these lines, either in full or upon any 
Phase. No information is available for many of the problems; however, 
with work being done at the present time on properties of metals at 
temperatures at and near the melting point, some hope is held that in 
the future information will be available concerning these properties. 

The process involves pouring molten steel (with a minimum amount 
of superheat to allow good fluidity) inte a vertical, water cooled, 
metal mold. At the start of the process the mold is plugged at the 
lower end until sufficient freezing takes place to allow adequate 
strength to permit pulling the plug and continuously processing the 
resulting billet. 


An analysis involves investigation of various phases which are 
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outlined below and discussed later in this paper: 


(a) 
(b) 
(c) 


(a) 
(e) 
(f) 
(g) 
(h) 


(4) 


Temperature distribution at any point in the ingot. 
Investigation of limiting water flow rates, and 
Prediction of accompanying coefficients of heat transfer 
between the ingot and the mold, 

Possibility of frictional drag exerted on the mold by 
the ingot. 

Maximum velocity of ingot travel. 

Load capacity of the ingot. 

Determination of ingot breaking point. 

Investigation of possible surface boiling on water side 
of the mold. 


Investigation of surface tension and surface film effects. 





HISTORY 


The greatest amount of work on continuous casting of steel has 
been done by Edward R. Williams and the Williams Engineering Company. 
His solutions to the problems have been of an experimental nature, 
and to the best knowledge of the writer, no attempt has been made 
toward an analytical investigation. 

The Williams process accomplishes the necessary heat removal by 
providing a very thin mold wall of high conductivity material (com- 
mercial drawn copper or brass tubing) through which heat can rapidly 
pass from the ingot, and by removing this heat extremely rapidly from 
the outside of the mold wall by means of water or other cooling fluid 
passed at high velocity along the surface. By removing the heat 
rapidly from the outer surface of the mold wall, a high temperature 
gradient is produced in the mold wall, which permits fast removal of 
heat from the ingot. It is claimed that this instantly freezes the 
molten metal into a skin which shrinks away from frictional resistance 
and is sufficiently strong to prevent fracturing as it is moved along 
the mold wall. 

The cooling water in the mold passes through passages which have 
a very small clearance, the velocity being quite high. It has been 
found that a relatively small amount of water is sufficient to satis- 
factorily remove the heat from a steel ingot to enable it to be poured 
and withdrawn at a speed of seven or more feet per minute. 


As the ingot emerges from the bottom it passes through a series 
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of sizing rolls that prevent bulging of the walls by ferrostatic pres- 
sure from the molten interior, Water sprays surround these sizing 
rolls and the ingot is quite cool when it gets down to the pinch rolls, 
Williams frequently uses mold lubrication to increase speed of cas- 
ting. The lubricant is introduced at the top along the mold wall, 

the quantity being controlled by the speed of the pinch rolls drawing 
the continuous billet out. 

The Williams machine has been given no steady trial over a long 
ran on steel; however, the experimental ingots produced have, for the 
most part, satisfactory structure. 

A number of makers of carbon steel who habitually handle small 
heats have shown interest in the process. It is also conceivable 
that those makers of high alloy steels and certain types of stainless 
steels may see in this process a possible means for eliminating present 
high scrap losses and a more flexible method of handling small tonnage 
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DETERMINATION OF THE TEMPERATURE AT ANY POINT IN THE INGOT 


In the following investigation certain assumptions are made to 
simplify the problem and bring it within the realm of possible solution. 
These assumptions are as follows: 

(a) Thermal conductivity does not vary with tem 
perature. 

(b) Wo recognition of the heat of fusion of steel 
is made in the temperature distribution in 
the ingot. 

(c) Convection in the molten metal is neglected. 

It is realized that the first two assumptions will cause appreciable 
error in the final results; however, unless these assumptions are made 
the problem would be beyond the scope of this paper. The assumption 
of the absence of convection in the molten steel is justified since it 
will become "mshy" upon entrance into the mold and resist any ap~ 
preciable turbulence. 

For purposes of this analysis a cylindrical ingot is used with 
all results in cylindrical coordinates. Consider first an elemental 


volume as follows; 
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The partial differcntial equation expressing the heat conduction 
for this problem is: (refer to Appendix "A" for derivation). 
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where co is the "thermal diffusivity" of the material, 


x= bbc, 


The thermal diffusivity is a property of the substance in which 
heat is transferred and stored. It is the ratio of a quantity which 
18 proportional to the heat conduction (numerator) and to a quantity 
which 4s proportional to the heat storage (denominator). Heat unites 
cancel and the dimensions become ft*/nr. 

Assuming it, is possible to express the solution of the above 
equation as the product of 3 functions, one of which is a function of 
fy alone, one a function of "f" alone, and the other a function of 


"3" alone, then 


te Re) BQ) Ze) 
The solution (as shown Appendix *BY) 4s 


t= 3.69 |e Veh 


where "a" is 2 constant 


C 


oo a 


Note: Sinee the temperature distribution is a function of 
"yx" and "s" only, terms including "() drop out. 


6. 
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The ingot is eonsidered to be ea semi-infinite cylinder with the 
following boundary conditions: 

az. t= §©) 220, OL KR 

BR tho , zar°, Fe 

3. t= F@ , o<2<co, c= 

From the second boundary condition it is seen that for too 


atzenr, Cs must be zero. The equation then reduces to 


se AF Gs) an (= a VG) +0") 


Solution for the above is made by assuming that the temperature "t" is 
the sum of two solutions, that is 
#£2t, +t, where 1256) 220 OSes 
t, 2%, (2,r) such that when ih O osaeo, «= R 
t= t, (zr) guch that Posie OstsK 
t,7 Fle), oszseo, F=K 
Solutions for t, and t, are taken as modified forms of equations 
31 and 32, pp 192, Carslaw and Jaeger, Conduction of Heat in Solids or 
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where 
N a 
a Ni 
Ky ax Ge) nth 
and the a, are the positive roots of J(Ra,) = 0. 


The total solution then becomes the sum of equations (1) and (2) above. 
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Reducing the latter part of equation (2) to more usable form, we have, 
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The temperature distribution in the he thus becomes 
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For actual temperatures the functions f(r) and FA mist be 
sapplied from physical aspects of the problem. Consider, for steady 
flow, the following assumptions: 

f(r) 2 Tp (pouring temperature) 
K)s ne “where L is the length of the mold. 
Substituting these values and integrating we have 
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Since K,. is negative for all values of a,, thie will converge very 
rapidly. 
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INVESTIGATION OF MAXIMUM COOLING WATER FLOW RATES, RELATION TO HEAT 
TRANSFER COEFFICIENT, AND EXPRESSION FOH MOLD LENGTH 


For this phase of the problem consider the magnitude of the 
heat transferred from the ingot to the cooling water with the as- 
sumption that the velocity of the ingot is such as to cause final 
freezing at the instant the ingot emerges from the mold. Since the 
mass rate of flow of the metal is constant, the heat transferred in 


unit time for a mold of length L is equal to 


Q2m(hy, - hoy) Btu/sec 


. war 
where m= sar Vy Vv, = velocity, ft/nr 
D = diameter ingot 
& = density steel 


fT, pouring temperatur 
H = heat of fusion 
h Cc_F Cesz Specific heat 
warhdg Tt « mean temperature of 
ingot cross-section 
at mold exit 


Substituting these values into the above equation, the total heat 


transferred is 


Q= Ts | 4 +O 7)| Blu |e, 


Dividing through by T(DL, we have 
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Since "q" and "V," are independent there is obtained an expression 


for the length. of the mold in terms of the heat transferred. 


lat +¢, Gr -7)| 


From the basic law of heat conduction for a hollow cylinder, 
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the heat transferred from the ingot to the cooling water is 
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where 


rj = radius of ingot 

ro = radius of mold, water side 

t] = average surface temperature of ingot within 
mold 

tg s average water temperature in the mold 


The term "U" is a composite or "overall" coefficient of heat transfer 


and is defined as followsi 





Pie \ 
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te ie ee on 


wh ere 
hy = — coefficient of heat transfer, water 
si 
k, = thermal conductivity of mold 
hy s sarface coefficient of heat transfer, ingot 
to mold 
1 s mold wall thickness 


from this, it is seen that 
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Substituting this into the equation of heat conduction and solving 


for h; 5 we have 
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For liquids with Reynold's numbers exceeding 10,000 the expression 


for hy is given by the Colburn equation 


(equation hd, pp 168, McAdams, Heat Transmission) 
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Applying this data to a specific problem, consider the fole 


lowing data: 


Diameter of ingot 
Pouring temperature, T 
Length of mold, L 

ew coning temperature 


Heat of fusion of steel 
Specific heat steel 
Mold thickness 

ky/2 


l, inches 
2800°F 

B Pt. 

27600F 

256007 

65 Btu/# 
0.1075 Btu/# 
0.5 inch 

1650 B/hr.ft.F 
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From this 


211.1 x10? Btu/hr 
= 3654 x 109 Btu/nr, ft? 


Using this value of "q" in the equation for the determination of 

the surface coefficient of heat transfer (ingot to mold), h, can be 
determined. With this value and values for h, and ky/l, an "overall" 
coefficient for heat transfer for varying Reynold's numbers is 
calevlated. Curves of hy, hy, k,/l and U versus Reynolds numbers 
are shown plate I, pp 14. 

The results of these curves show that the overall coef- 
ficlent of heat transfer is dependent upon the surface coefficient 
between the ingot and the mold, and for this particular set of con- 
ditions, water flow rates in excess of Heynolds number 5 x 104 would 
be impractical. Such determinations can be made for varying sets 
of conditions giving a family of similar curves, and from these flow 
rates, mold length, and surface coefficients of heat transfer can 
be predicted. These results are not conclusive as much depends 
upon reasonably accurate determination of temperatures within the 
ingot. It dees, however, show that maximum cooling water flow rates 
are dependent upon the amount of heat that can be transferred from 
the ingot to the mold, which in tum is dependent upon the correspond- 


ing surface coefficient. 
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CONCLUSIONS 


1. As pointed out previously the solution for the temperature 
distribution is only approximate because the heat of fusion of the 
metal has been neglected in the derivation. The solution, however, 
is useful in heat transfer problems for determining temperatures in 
cylindrical bodies moving with constant velocity, and also in 
predicting temperature distributions in contimous casting problems. 
Solutions for different velocities can be modified by proper selection 
ef the functions f(r) and Fe) which can be fairly accurately pre- 
dicted from experimental data. 


2. In the investigation of maximm water flow rates, and their 
relation to surface coefficients of heat transfer between ingot and 
mold, it is shown that the over-all coefficient of heat transfer is, 
in the limit, a function of the surface coefficient between the billet 
and the mold. From the basic law of heat conduction and for a given 
sise ingot, it is concluded that to increase the heat transferred some 
means must be provided to increase this surface coefficient. In the 
Williams process mold lubrication is used, which, in addition to re- 
ducing friction a resistance, increases the surface coefficient of 
heat transfer by making the heat transfer across the space more nearly 
that of pure conduction. It is recommended that for the lubricant 
an oil of high thermal conductivity be used - the optimum between 
the lubricating effect and the increased heat transfer coefficient is 


left for future investigation. 
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After determination of the maximam heat transfer rate, the 
optimum length of the mold can be determined from the relation 


It is to be noted here, however, that the above equation is based on 
the assumption that final freezing takes place just as the ingot 
emerges from the mold. More realistically, this final freezing point 
takes place some distance below the mold. That is, the mold length 
need be only sufficient to establish a "skin" around the molten metal 
which is strong enough to prevent rupture either from excessive load 
or from ferrostatic pressure from the molten interior. Investigation 
of minimum mold length is recommended for further work. 
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RECOMMENDATIONS FOR FURTHER WORK 


Since the technique of continuous casting is complex, and 
since time limitations precluded analysis of all phases, the following 
recommendations are made for considerations 

1. Frictional drag exerted on the mold by the ingot 

Although it is claimed that the molten metal instantaneously 
freezes upon entry into the mold and shrinks away from frictional 
resistance (The Iron Age, April and April 11, 1940), it atin highly 
improbable. One of the laws of sliding friction is that the friction 
is proportional to the force pressing the surfaces together. Thie 
force is increased by the ferrostatic pressure from the molten 
interior, which, it is believed, will overcome the shrinking due 
to freezing of the metal. Adding a lubricant to decrease the coef- 
ficient of friction will increase the maximum load capacity and hence 
the maximum allowable velocity. 

Investigation of the mechanical friction is simplified be- 
cause, within limits, sliding friction is independent of the area of 
contact, and (except at start) of the speed of relative motion. 

The problem, therefore, resolves itself to a balance of forces be- 
tween the ferrostatic pressnre and the shrinking due to cooling. 
2. Maximim Velocity of ingot travel. 

The velocity of the ingot varies directly as the amount of 
heat which can be satisfactorily removed from the steel. With suf- 
ficient heat removal from the mold a skin will be developed such that 
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the ingot can withstand the load imposed by the pinch rolls. The 
maximum velocity will be determined by the maximum heat removal 
consistent with strength characteristics of the skin of the billet. 


3. Load capacity of ingot and determination of 
ingot breaking point 


The load capacity of the ingot is a function of the load 
capacity of the skin formed in the mold. The load will be caused 
by surface tension and surface film effects of the steel and friction- 
al resistance to the passage of the ingot. This load capacity will 
reach its maximum when such velocity is reached that the maximum 
heat removal will be below that necessary to develop sufficient skin 
thickness and strength. Since the ferrostatic pressure increases 
down through the mold, the frictional resistance will also increase, 
and hence some point will be reached where the load exceeds the 
ultimate strength of the billet. Experimental evidence has shown 
that the breaking point will occur some distance down the mold, but 
the exact location has not been determined. 


. Investigation of possible surface boiling on the water 
side of the mold 


With very high rates of heat transfer and high temperature 
gradients produced in the mold wall there are excellent possibilities 
ef surface boiling in the cooling water. With high feed rates, a 
small fraction is vaporized, and the boiling-side coefficients may 
be roughly twice those predicted from straight water cooling. These 
increased coefficients may be due to increased turbulence in the film 


near the wall. Investigation along these lines can be very instru- 
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mental in increasing heat transfer rates. 

S. Surface Tension and Surface Film Effects 

A high value of surface tension has the effect of increasing 
tae pressure required to cause the metal to enter and flow along a 
narrow channels but the direct effect of the true surface tension 
is almost overshadowed by the effect of surface films, the influence 
of which is equivalent to that of a greatly increased surface tension. 
It has been found that surface oxide films affect the flowing power 
of metals with the apparent surface tension of some metals increased 
three timess Steel has a very high value ef surface tension. 
(Briggs, Metallurgy of Steel Castings) 

Both effects will add te the resistance to flow and decrease 
the maximum load capacity of the ingot. Although it is thought to 
be considerably small compared to frictional resistance, further 
investigation should be initiated to determine its effects. 
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APPENDIX A 
Equation of Conduction for a Moving Cylinder 





Heat conducted through "r* face is 
wt dade - | Fees - dQ, 


Heat conducted through the "r+ dr" face is 


kdods © rdt + kdddes” Cdr dy tkddde = Se. JQ, 
Heat conducted through "&* face 
Ini NR 7 = - ¥-400,, 
r@ 


Note: Temperature gradient in Pdirection is the 
change of temperature with respect to the 
distance represented by the angled? which 
distance also depends on the radius r. 


Heat conducted aes (Pp + df " face 
er Sear CAE elves) Lay *@ de 


el. 
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sh) be» iodil ise: = ath 4 ne ababa 












» st 
st ore 
eiear. meee, "3. at 
Since sea)barg | * " Hea) ° ved ae 
v 
then | IO,9 = 
i Ny dat 9) 2 aE 
kdrde fg) torkdedadd Sa dt 20 
Heat conducted through "2" face 
k= £ rdadede > dO. 
Heat conducted through "z +dz"* face 
k 2 rdddedt + kdddr & or 2 de | de 
eee SE rife dt ipreerie * Fen oe 
Heat gained by the elemental volume due to conduction is equal to 
7d oe " dQ 
which equals 


br eidodedadr 4 \ dards de + : 3 alaldedt mee t de dade 


Heat stored in the differential body 
da.: e Ce (cdadzedr dt 


However the temperature increment of a particle of the fluid in 


time interval dT becomes 


Lat “ 
aie x de + ry dO 7 a 582 + St dr 
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Gividing turough by dT 
at Ly bt, tvp th + yet - St 
AT dr z sa ST 
Sinee we are considering velocity in the sedirection and a condition 
of steady state, this reduces to 


dt= Vv, 28 de 


From this, then, the Heat stored in the body becomes 


da,= ie Ce (cdddede) eeu 


Since no heat energy is being developed in the body, the law of the 
conservation of energy reduces to 
Eeat stored = heat in — heat out 


= dQ, —- dQ, 


Substituting and dividing through ty rAOdrdzedT 2 
ot Xt s7é “ 
Ce Ve VS resales tos ye PYG din: 
introducing - term "thermal diffusivity” 


oy = eS 
(The thermal diffusivity is a property of the substance in which 
heat is transferred and ctored). 
Substituting this value, the equation becomes 
Ae Poses vt eat We ain 8 
st Y on rah = so” X 92 
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APPENDIX 3B 


SOLUTION FOR THE EQUATION OF CONDUCTION = TEMPERATURE DISTRIBUTION IN 


INGOT 


Assume it is possible to express the solution of the dif- 


ferential equation (Appendix A) as the product of three functions, 


one of which is a function of "r® alone, one a function of "¢* 


alone, and the other a function of "Z" alone, then 


t= Rr) BO Zw 


The equation takes the form, 


ZERYLZSR HRZSR + ZOSR- 


Vz re 
a ee, =O 


dividing through by RD-Z , transposing and equating to a constant, 


we have 


“dul 
ac ee 


1 
c *@ “= 
From the left side of the equation, write 


as + re ee _= 


which gives 


at 


e +p p=o 


| 
+> 
N 


the solution of which is 


Coie Camo i C,cos pO 


2h. 


co ae ) 


Q) 


(3) 


@) 





Equation (2) also gives 


re +#rkR R +(ra- ) R (Ss) 
let x = ar 
then since 

dR _ dR dx and dy ae 

ae.” dy dr dv 


equation (5) becomes 
CRark {x= rp) K=O (6) 
where = S@) = Sar) 


This equation is Bessel's differential equation (Hildebrands Advanced 
Calculus of Engineers,.pp 165-167), and the general solution for "p* 


an integer can be written 

R= To) +c, Kp) (7) 
where K, (x) is given as follows (Woods, Advanced Calculus, pp 288) 

Kp) = Spe) Into + ey 

42n 
ee @-nlxt Xx ei) y xt"? ie 

Where Py => Ly = = See \ Airs “nw 
Since equation (7) is applied to a solid cylinder, Bessel's function 
K, (x) becomes infinite at x = 0 (In 0. =-co), and it follows that 


Co = 0. Thus for the solid cylinder, the general solution is taken 


R= CI = CIpen (8) 
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Taking the remaining equation in our solution 


Z-4Z-dZ=0 @) 


Using the operator method of solution, the equation becomes 


[s-o-ee- 6 


The solution, therefore, of equation (9) becomes 


mt Mat 
£= C1 + Ce 


(10) 


(1) 


where y a 
= =o Feo 
a ) 
Vv Va" 
= 2 _ z t 
2 2, ) +o 


Hence for the cylinder the solution of the differential equation is 
given by substituting equations (lh), (8), and (11) into equation 


t= Re) Ba) Ze 


Con tS Tolar) |Csimpd +C, coop legre"| 


Since the temperature distribution is a function of "r® and "g® only, 


then, p = O, and the solution of the equation becomes 


t= Jar) loo oe] 
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